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We study the interaction of magnetic monopoles and domain walls in a model with SU(5)xZ2 
symmetry by numerically evolving the field equations. We find that the monopoles unwind and 
dissipate their magnetic energy on collision with domain walls within which the full SU(5) symmetry 
is restored. 
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The interactions of topological defects can have a pro- 
found effect on the outcome of phase transitions. The 
scaling of a network of domain walls and strings, and a 
distribution of magnetic monopoles, crucially depends on 
how the defects interact among themselves and with each 
other. Thus far attention has focussed on the interactions 
of walls with walls, strings with strings, and monopoles 
with monopoles. The cosmological importance of the in- 
teractions of walls and monopoles was highlighted in Ref. 
and it is this problem that we study in the present pa- 
per. 

Earlier work on the interaction of solitons and do- 
main walls (phase boundaries) has been carried out in 
the following contexts: (i) mutual interaction of domain 
walls 1^, (ii) Hc^ A-B phase boundaries and vortices 

, (iii) Skyrmions and domain walls [Q , and (iv) global 
monopoles and embedded domain walls in an 0(3) linear 
a model Here we will numerically study the interac- 
tion of gauged SU(5) monopoles with a Z2 domain wall. 
This is quite distinct from the earlier work since it looks 
at magnetic monopoles which necessarily include gauge 
fields. It is also the most relevant problem for the cos- 
mological consequences of Grand Unified theories 0] . 

The SU(5) model we consider is given by the La- 
grangian: 
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(1) 



where $ is an SU(5) adjoint scalar field, (a = 

1, ...,24) are the gauge field strengths and the covariant 
derivative is defined by: 



(2) 



and the group generators are normalized by Tr(TaTb) = 
5ab/2. The potential V{^) is the most general quartic 
potential but we exclude the cubic term in $ so as to 
obtain the extra Z2 symmetry under $ — > — $: 



F($) = -m^Tr*^ + /i(Tr$2)2 + ATr$* . 



(3) 



The parameters of the potential are chosen so that 
($) = 77diag(2,-3,2,2,-3)/(2yT5) with 77 = m/VV 
and A' — h + 7X/30. With this vacuum expectation 



value, the SU(5) symmetry is spontaneously broken to 
SU(3)xSU(2)xU(l). The desired constraints on the pa- 
rameters are: A, A' > 0. 

The magnetic monopoles in this model were discussed 
by Dokos and Tomaras Q except that also included the 
effects of a scalar field in the fundamental representation 
of SU(5). Here we do not have such a field. Yet the ba- 
sic construction of |^ goes through and the fundamental 
monopole is essentially an SU(2) monopole embedded in 
the full theory. The monopole solution has the following 
form: 



(4) 



where the subscript M denotes the monopole field con- 
figuration, 

T'^ = idiag(a^ 0, 0, 0) , = ^(0, 0, 1, 1, -2) , 
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(-3,-3,2,2,2) 
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cr" being the Pauli spin matrices, 

$° = P(r)x" , = M{r) , $^ = N{r) , (5) 



where r = \J + y"^ + z"^ is the spherical radial coordi- 
nate. The ansatz for the gauge fields for the monopole 
is: 



=^\^{^-K{r)) , (a =1,2, 3) , 



W,^ = 0, , (M 1,2,3). 



(6) 



In the case when the potential vanishes (the 
Bogomolnyi-Prasad-Sommerfield (BPS) case |@,||), the 
exact solution is known fOll: 
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er^ tanh(Cr) 



1) , K{r) 
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In the non-BPS case, the profile functions P{r), K{r), 
M(r) and N{r) need to be found numerically. We find 
them by using a relaxation procedure with the BPS so- 
lution serving as the initial guess. 

Depending on the parameters in the potential, it is 
possible to have different stable domain wall solutions. 
The domain wall across which $ ^ — $ is stable provided 
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At the center of this wall, <i> must necessarily vanish and 
so the full SU(5) symmetry is restored at the center of 
this wall. Certain components of $ do not vanish at the 
center of the domain wall solutions in this model for other 
values of parameters. In these walls, only a subgroup of 
the full SU(5) symmetry is restored in the center. We 
will only study the interaction of monopoles with walls in 
which $ = at the center in this paper. The interactions 
of other types of walls and monopoles will be discussed 
separately. 

The solution for the domain wall located in the xy- 
plane is 



V 



2^/T5 



tanh(crz)(2,-3,2,2,-3) , 



(10) 



where a = r]y/X'/2. 

When the monopole and the domain wall are very far 
from each other, the joint field configuration is given by 
the product ansatz: 



$ = tanh(70-(z - zq))^ 



M 



(11) 



where v is the velocity of the domain wall in the negative 
z-direction, 7 = 1/Vl — v'^ is the Lorentz factor and zq is 
the position of the wall. Here <I>j\/ denotes the monopole 
solution in eq. (^. The gauge fields are unaffected by 
the presence of the wall and are still given by eq. (||). 
In addition, the time derivative of the scalar field is also 
given by the product ansatze: 



$ = jav secli^(7(T(z — za))^M 



(12) 



Eqs. ( |ll[ ) and (g_l|) specify the initial {t — 0) conditions 
for the scalar field for a wall approaching a monopole 
with velocity v. The initial scalar and gauge field profile 
functions P, M, N and K (in the non-BPS case) are 
found by numerical relaxation. The field dynamics is 
described by the equations of motion following from the 
Lagrangian in (|l|). At first sight, there are 24 components 
of $ and 96 components of the gauge fields that need 
to be evolved. However, it is not hard to check that all 
the dynamics occurs in an SU(2) subgroup of the original 
SU(5). This then reduces the dynamical fields to a triplet 
of SU(2) and two other fields (i.e. a total of 5 scalar 
fields) and 3x4=12 gauge field components. Choosing 



the temporal gauge {Wq = 0) reduces the number of 
gauge field components to 9. 

Further reduction of the problem occurs since the ini- 
tial conditions are axially symmetric and the evolution 
equations preserve this symmetry. The angular depen- 
dence in cylindrical coordinates can easily be imposed 
on the scalar field. For the gauge fields it can be ex- 
tracted by using the fact that the covariant derivatives 
of the scalar field must vanish at large distances from the 
monopole. This then leads to the following ansatz for the 
5 scalar and 9 gauge fields: 

$1 = /i x , $2 = /i y , ^j. = h z 
2 I 
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where the fi (i = 1,...,8) are functions only of i, p = 
x"^ + y'^ and z. We have explicitly checked that this 
ansatz is preserved by the evolution equations. So now 
the problem is reduced to one in 8 real functions of time 
and two spatial coordinates. 

An attempt to numerically solve the 8 equations of mo- 
tion directly in cylindrical coordinates failed due to nu- 
merical instabilities that developed within the time scale 
of the simulation. An analysis showed that the problem 
was due to large numerical errors in evaluating the deriva- 
tives in cylindrical coordinates. This shortcoming of us- 
ing cylindrical (and spherical) coordinates in numerical 
work is well-recognized and the authors of [|l0| have pro- 
posed a solution that we have successfully implemented. 
The idea is to solve the problem, not in two spatial di- 
mensions like the pz-plane, but to solve it in a thin three 
dimensional slab whose central slice is taken to lie in the 
xz-plane and with only 3 lattice spacings along the y di- 
rection. Then Cartesian coordinates can be used to solve 
the equations of motion in the y = plane, thus mini- 
mizing numerical errors. On the y ^ lattice sites the 
fields are evaluated by using the axial symmetry of the 
problem. This scheme improved the numerical stability 
of our staggered leapfrog code dramatically and allowed 
us to observe the monopole and wall for a sufficiently 
long duration without the development of numerical in- 
stabilities. 

We have evolved the initial wall and monopole config- 
uration with several velocities. The numerical results of 
the simulation with 1; = 0.8, h = — A/5, A = 0.5 and 77 = 1 
are given in the figures and clearly show that the energy 
of the monopole dissipates after the passage of the wall 
(m = ri\/X' and e can be scaled out of the problem). The 
final snapshot shows that the energy in the scalar field is 
located entirely on the wall and the magnetic energy is 
along and behind the wall. 

When the domain wall moves close to the monopole 
the latter is pulled toward the wall. This signals the 
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presence of an attractive force between the two defects. 
Such a force is expected from energy considerations ||] 
and has been observed in the 0(3) linear a model studied 
in|]. 

We have estimated the time it takes for the monopole 
to be destroyed as a function of different wall velocities. 
The topological winding density of the monopole is a 
delta function in space, while the topological winding 
is a discrete number and so these quantities cannot be 
used to estimate the dissolution time. Instead, the mag- 
netic energy density provides us with a continuously vary- 
ing quantity that we can track in the simulations. (The 
scalar energy is not suitable since both the wall and the 
monopole contribute.) We have chosen a cylindrical vol- 
ume around the monopole with axis along the z-axis and 
computed the magnetic energy within this volume as a 
function of time. This gives us the rate at which the mag- 
netic energy escapes the cylindrical volume surrounding 
the monopole. There are two stages in the dissolution 
process - first the monopole gets absorbed by the wall 
and then the magnetic energy starts spreading in the di- 
rection along the wall. The time taken for the monopole 
to be absorbed by the wall is measured from when the 
cores overlap (corresponding to a sharp drop in the mag- 
netic energy and a sharp rise in the electric energy) to 
when the magnetic energy starts moving together with 
the wall in the z-direction but spreading along the wall. 
In our simulations, we observe that the absorption time 
is approximately equal to the width of the wall in the 
monopole rest frame, that is tabs ^ (70")^^ • During the 
second stage the remaining magnetic energy is confined 
to the wall and is escaping the cylinder as it travels along 
the wall. In the rest frame of the wall the duration of the 
second step is independent of the wall velocity but takes 
longer in the initial rest frame of the monopole by a factor 
that is well- accounted for by time dilation. 

At high wall velocities, Lorentz contraction results 
in the wall appearing to be much "thinner" than the 
monopole. We have tried velocities as high as w = 0.99 
and have not observed any qualitative alterations from 
the dynamics at lower velocites except for the obvious 
time dilation of the propagation of the excitaions along 
the wall. 

The above results have not changed as we varied the 
parameters in the potential within the range in which a 
$ — $ domain wall solution is stable. To better under- 
stand this independence of the dynamics on the param- 
eters let us examine the parameter space itself and the 
characteristic length scales that are involved. The stabil- 
ity of the domain wall is determined by the ratio h/\. 
The value of rj sets the fundamental time and length 
scales of the simulation and can be used to adjust the 
lattice grid spacing to the "widths" of the defects. There 
are three relevant length scales in the problem: the sizes 
of the scalar and vector cores of the monopole, rs and 
ry, and the thickness of the domain wall, r^w. These 



are the sizes of the regions in which the fields deviate 
significantly from their asymptotic values. We can es- 
timate r5, ry and row numerically by finding the dis- 
tance at which the corresponding fields become an expo- 
nent closer to their respective asymptotic values. From 
dimensional agruments rs and ry are approximately in- 
versely proportional to the masses of the scalar and vec- 
tor bosons making up the monopole, ms = Tj^/X and 
mv = ■\/5/12 erj. The thickness of the domain wall 
is given by a^^ — The interaction between 

the monopole and the wall depends on the relative val- 
ues of rs and row- The dimensional arguments give 
rs ^ row while numerically we find rs < row for 
all values of h/X that lead to stable domain walls (see 
eq. (|)) and A £ [0.01,100]. For sufliciently high val- 
ues of A it is possible to have the vector width of the 
monopole to be larger than the domain wall width, i.e. 
ry > rjjw > rs- Our simulations show that the walls 
still sweep the monopoles. We will study the many dif- 
ferent kinds of domain walls that can occur in the model 
and their interactions with monopoles separately. 




FIG. 1. The first panel shows the potential energy den- 
sity in the xz-plane for the magnetic monopole and domain 
wall where h = —A/5, A = 0.5 and the wall velocity is 0.8c. 
The second panel shows the corresponding magnetic energy 
density (proportional to B°'1). 
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FIG. 2. As in Fig. hi at an intermediate time step. 



As is clear from Fig. ^, the final state of the wall is 
different from that of the initial wall. The destruction of 
the monopole has left a residue of scalar and magnetic ex- 
citations on the domain wall that are propagating along 
and behind the wall. 

The dissolution of magnetic monopoles by domain 
walls implies that the number density of magnetic 
monopoles will fall off faster than if there were no do- 
main walls. The cosmology of such a system of walls 
and monopoles has been discussed in Q where it was 
argued that such interactions might resolve the cosmo- 
logical monopole over-abundance problem. Similar in- 
teractions between strings and domain walls would affect 
the cosmological implications of cosmic strings. The nu- 
merical techniques presented here can also be used to 
study the interactions of walls and (global) monopoles or 
vortices in other systems. 
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FIG. 3. As in Fig. m at the final time step. 
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